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^ ■ Recently, Kovner and Lublinsky proposed new small-x QCD evolution equations valid when the 

. gluon density inside the target is low. The key element of their construction is the Wess~Zumino 

term which ensures the non-commutativity of valence charges. In this paper we clarify the origin 
and significance of this term by showing that it can be naturally incorporated in the effective 
. theory of Color Glass Condensate. We also reexamine the renormalization group description in 

\ the high density (JIMWLK) regime. 
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1 Introduction 



in 
o 

Over the past few years, considerable effort has been made [H3] towards understanding 
^ . the high energy QCD evolution equation beyond the Balitsky-Kovchegov (BK) equation 

[14,15], or the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) 
\ equation [16,17]. Interest in the problem was sparked by several seemingly unrelated 

;h ' observations [1-4, 18] regarding the insufficiency of the BK-JIMWLK equation, but they 

all point to a necessity to include new types of Feynman diagrams called Pomeron loop 
diagrams. 



Fig. 1(a) is a typical diagram summed by the JIMWLK equation. The upper blob repre- 
sents an energetic hadron (the "target" ) moving in the positive z direction. At very high 
energy, the target behaves like a weakly coupled many body system of small-x gluons 
commonly dubbed the Color Glass Condensate (CGC). In one step of quantum evolution 
towards smaller x, arbitrarily many t-channel gluons inside the target recombine into two 
gluons. This is the gluon saturation phenomenon [19-22] which plays an important role 
for the unitarization of the BFKL Pomeron [23]. 

Attention is currently focused on diagrams with arbitrarily many gluon legs as represented 
by Fig 1(b). Viewed from the target side, these diagrams describe the gluon splitting or 
the gluon number fluctuation inside the target [4]. They are particularly important when 
the target is dilute, and are responsible for the eventual formation of a dense system. On 
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(a) (b) (c) 

Fig. 1. Quantum evolution of the target (represented by a blob): (a) Gluon recombination (b) 
Gluon splitting (c) Both the gluon recombination and the splitting 

the other hand, the same diagrams, when viewed from the projectile side, clearly describe 
the gluon recombination process in the projectile. In either interpretation. Fig. 1(b) is 
the leading diagram in the kinematic regime which is complementary to that considered 
in the JIMWLK equation. A set of evolution equations including the vertex of Fig. 1(b) 
with four gluon legs below the rung has been proposed as the first step beyond the 
JIMWLK equation [4,5, 12]. [See also earlier works [24].] These results were obtained in 
the frameworks of Mueller's dipole model [25] and its alternative formulation as a color 
glass [26]. In this model, the relevant gluon sphtting process (or better be called the 
Pomeron splitting process in the large-iVc approximation) is naturally described by the 
dipole splitting and subsequent gluon emissions. Combining the dipole model with the 
JIWMLK formalism in the dipole sector [27], one can investigate more general problems 
which involve both the Pomeron splitting and the Pomeron recombination (Fig. 1(c)), 
i.e., the Pomeron loops [9,24]. 

Meanwhile, in a very interesting paper [6] Kovner and Lublinsky have derived evolution 
equations in the dilute regime to all orders in the gluon legs without using the large- 
Nc approximation. The kernel of the evolution equation was shown to be dual to the 
JIMWLK kernel [7, 10, 11] essentially reflecting the symmetry between Figs. 1(a) and 1(b). 
However, identification of states ( "observables" ) on which this kernel acts does not follow 
from simple duality considerations. The nature of the evolution equation in the dilute 
regime is peculiar; one has to treat color charges as non-commutative operators. ^ Kovner 
and Lublinsky overcame this problem by introducing a Wess-Zumino term and an extra 
coordinate ("ordering variable"), after which the charges can be treated as commutative. 
The emergence of this term at first sight is somewhat mysterious, as it is rarely discussed 
in the context of high energy QCD. Still, the fact that the duality and the relationship 
between the wave function approach [6] and the effective action/Hamiltonian approaches 
[10, 11, 21] become manifest only in this "commutative world" calls for a formulation which 
explicitly includes the Wess-Zumino term. [See, however, [8].] 

In this paper we fully clarify the origin and significance of the Wess-Zumino term in the 
high density (JIMWLK) regime. This may be unexpected as the issue of non-commutativity 

^ As shown in [12], this problem disappears in the dipole model in the large- A^ approximation, 
and one can check the consistency with the previous results in [4]. 
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seems least important in this regime. On the contrary, we shall show that the Wess-Zumino 
term is quite useful and can be naturally incorporated in the JIMWLK formalism as a 
part of the source term. This observation has led us to more detailed treatment of the 
renormalization group description than in the literature [16, 17] in the sense that we 
discuss renormalization of the source term. It should be emphasized that, although in 
this paper we focus our attention to the JIMWLK regime, our approach is general and 
therefore applicable to both the dilute regime (Fig. 1(b)) and the Pomeron loop regime 
(Fig. 1(c)) where the importance of the Wess-Zumino term has been first recognized. 
Indeed, an essential difference between effective theories in these regimes and the current 
JIMWLK formahsm [16, 17] is that one has to introduce the light-cone time coordinate 
x"'"(=ordering variable) in the former [11]. By reformulating the latter with a; ''"-dependent 
classical charges, we can address the high energy evolution in different regimes in a single 
framework. In this context the Wess-Zumino term necessarily arises as a consequence of 
gauge invariance. 

In Section 2, after a brief review of the JIMWLK formalism we propose a new source term 
which explicitly includes a Wess-Zumino term. We show that our source term correctly 
reproduces the induced charge upon quantum evolution. In Section 3, we perform the 
renormalization group evolution directly at the level of the JIMWLK functional integral. 
We conclude with Section 4. 



2 JIMWLK formalism with the Wess-Zumino term 



2.1 The Color Glass Condensate 

The color glass condensate (CGC) formahsm [16, 17, 20, 22] is an effective theory of gluon 
saturation at a small Bjorken parameter x = A~^/P~^, where P"*" and A"*" are the light- 
cone momenta of the right-moving parent hadron (the target) and the small-a; gluon of 
interest, respectively. Partons with momentum fraction larger than A+ are described as 
the external source p, while the gluons with momenta A+ are described as the classical 
field created by the source according to the Yang-Mills equation 



We use the notation x = {x~,x^) where a;* {i = 1,2) is the two-dimensional transverse 
coordinate. In the saturated regime (JIMWLK regime), p is effectively static (i.e., 
independent) and parametrically of order ~ 0{l/g). In the hght cone gauge — 0, the 
solution to Eq. (2.1) is given by 



(2.1) 



Af" = S^^'B', 



= -Ud'U\ 
9 



(2.2) 
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where C/ is a Wilson line in the x direction 

U\x) = Pexp (^g j"" dz-a{z-,x')^ . (2.3) 

a is the classical solution in the Coulomb gauge 

= (2.4) 

and is related to the Coulomb gauge charge p = WpU via the plus component of the 
Yang-Mills equation 

-Via = p. (2.5) 

Observables O (e.g., scattering amphtudes) are first calculated for a given background 
field Eq. (2.2), 0[p\, and then averaged over p with the weight function Wr[p] 

{0)r = J DpWr[p]0[p], (2.6) 

where r is the rapidity r = In 1/x = In P+/A+. Wr satisfies the JIMWLK equation [16, 17] 
which is a renormalization group equation in rapidity 



d 

-K-^Ap] — — -f^JIMWLK 
OT 



s 

a, — 



WM- (2-7) 



The precise form of the JIMWLK Hamiltonian i^jiMWLK is irrelevant here. [See, however, 
Eq. (3.34).] Suffice it to say that it is quadratic in the functional derivative 5/5p (cor- 
responding to the two gluon legs in Fig. 1(a)) and all orders in a (corresponding to the 
merging gluons above the rung in Fig. 1(a)) in the form of the Wilson fine Eq. (2.3). 

The derivation of Eq. (2.7) proceeds with the following steps [16,17]: (i) Start with the 
QCD functional integral in the light-cone gauge 

j DpWM DA''5(A+)e^^YM[^''l+^5w[^^pl^ (2.8) 

where the subscript r in the integral means that the gauge fields contain only modes 
with j9+ < A+ = e^^P+. Modes with p+ > have been already integrated out. S'w is 
the source term which gives the right hand side of the Yang-Mills equation Eq. (2.1). (ii) 
Expand the gauge field around the solution to the Yang-Mills equation 

^ ^P-iQi + + (2.9) 

where is the semihard field with the momentum fraction A+ > p"*" > 6A+ {h <^1) and 
5A^ is the soft field with p"*" < feA"*". (iii) Functionally integrate out the semihard field 
a^. (iv) Show that, in the leading logarithmic approximation (LLA), i.e., keeping terms 
proportional to In 1/6 = UgSr (a^ = g^/in), the effect of the integration is absorbed 
by the renormalization of the weight function Wrlp] — > Wt+st[p + ^p], where Sp is the 
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additional charge induced by the semihard field. In this paper we perform this program 
in the most direct way, including the renormalization of the source term 5'w- 

For this purpose, first we must know (S'w- The original definition in [16] is 

^w=^/ dftr[p(f)ty(f)], (2.10) 

where 

W{x) = Pexp (^ig p dx-^A-{x+,x)T''^ , (2.11) 

with T°- the color matrices in the adjoint representation. [Hereafter W and W denote 
Wilson lines in the fundamental and adjoint representation, respectively.] Eq. (2.10) is 
invariant under x"*" -independent gauge transformations. Upon quantum evolution, one 
shifts A~ ^ a~ + 5A~ [Note that = for the classical solution.] and expands the 
exponential. The coefficient of —SA'^ defines the induced charge 6p[a~] to be included in 
the classical theory at rapidity t + 6t = lnP+/6A+. Several different choices of can 
be found in the literature [17, 32]. The non-uniqueness of S'w simply means that they are 
effective actions. Although they all lead to the renormalization group equation Eq. (2.7), 
it is not clear whether intermediate calculations go hand-in-hand with the corresponding 
perturbative QCD calculation. Moreover, when we consider the all order effect of A~ in 
the dilute regime, the choice of S'w is crucial to obtain correct evolution equations. 

Below we construct another Sw that is most suited for our purpose. This is not in the 
form of an effective action, but contains its own dynamics due to the presence of the Wcss- 
Zumino term. For the sake of simplicity and clarity, throughout this paper we consider the 
color SU(2) gauge group. ^ We describe the target hadron as a bunch of energetic "valence 
partons" having longitudinal momenta p"*" > e~'^A+. They carry various representations 
of SU(2) {J — |,1,|,---) at varioTis transverse coordinates x\ and are also distributed in 
a narrow strip in the x~ direction with the width |Aa;^| ~ 1/A+ (Fig. 2). By a deliberate 
choice of the light-cone prescription, one can restrict the support of the charges to the 
region x~ > [17]. 

Let us focus on a single quark at the coordinate x inside the target. In the eikonal ap- 
proximation, the propagation of this quark is described by the amplitude 

Z^{4^aW''\x)^l^,), (2.12) 

where W is the Wilson line Eq. (2.11) in the fundamental representation T" r"/2 (r" 
a = 1,2,3 are the Pauli matrices) and ijj is the Dirac spinor of the quark. With the use of 
a formula for an open Wilson line W"'^ derived by Diakonov and Petrov [29], Eq. (2.12) 
takes the form 



^ The form of the Wess-Zumino term is a bit simpler for SU(2) than for SU(3). The latter can 
be found in the literature. See, for example, Ref . [29] . 
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Fig. 2. Eikonal propagation of the valence charges. The charges have support only at positive 
x" . 



J DS±ooW[Soo,S.oo]J^°^ DS{x^) J DA^ 
X exp (^iSyM + igJ J 

J Jp-oo J 

X exp [iSym - i J dx+p{x+)A-{x+) + iSwz[p{x^)]^ , (2.13) 

where S{x~^) is a SU(2) matrix {S±qo = S{x^ = ±00)), J = 1/2 for the quark fundamental 
representation. S'wz is the Wess-Zumino term (also called a geometric phase, Berry's 
phase, Polyakov's spin factor) [28, 29] 

/oo 
dx+iY[T^Sd-S\ (2.14) 
-00 

where d~ — d+ = d/dx'^. In the second equality in Eq. (2.13), we switched to the p- 
representation defined by {A~ — A~ ^) 

p«(x+) = -^tr [t^St'^S^]. (2.15) 

The constraint p"p" = g^J"^ is implicit in the measure Dp. For a Wilson /oop, S'wz can be 
written as a functional of p by introducing yet another coordinate -u (1 > -u > 0) (which is 
the radial coordinate of the disc spanned by the loop) and extrapolating p{x^^ — > p{x^ ^ u) 
such that p(a;+,n = 1) = p(x+). But for an open Wilson line, this is not possible in 
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general.^ Nevertheless, we use the notation S'wz[p] with p and S related by Eq. (2.15). 
The "weight function" W^[pooiP-oo] (not to be confused with the Wilson line W) contains 
Wigner rotation matrices [29] and quark spinors. It is in general complex and depends on 
S (or p) only at x'^ — ±00. Note that the sum of the last two terms in the exponential of 
Eq. (2.13) 

gJii[T:iSA-S^] + iJtr [r^Sd-S^] (2.16) 
is invariant under the following gauge transformation 

A- = U^A'U + -U^d-U, S = SU (2.17) 

Under this transformation, p transforms as 

p = U^pU. (2.18) 



In fact, Eq. (2.13) is a path integral formula for a spin [28]. The exponential factor 

H{x^) = -gA-{x+)^, (2.19) 

of the fundamental Wilson line in Eq. (2.12) may be regarded as a Hamiltonian for a 
non-relativistic "spin" r" immersed in a time-dependent "magnetic field" A~{x'^). The 
Wess-Zumino term is nothing but the kinetic term (the pq term in Eq. (2.20)) which 
arises when one goes from the Hamiltonian to the Lagrangian^ 

Pexp (^-i J dx+H^ ~ J DpDqexp (^i J dx+{pq - H)^ . (2.20) 

Generalization to an arbitrary set of valence partons (with longitudinal momenta p"*" > 
e~'^A"'") can be done simply by replacing 

J dxJ(x), (2.21) 

and endowing Vrr[p±oo(3^)] with the information about the distribution of partons and 
their group representations. This can be done along the lines suggested in [33]. In this 
way we are led to the following source term {x = (x+, x)) 

Sw[A-, P] = - / dxp{x)A-{x) + Swz[p{x)i (2.22) 



J 



^ Wc could consider an overall color-singlct target and make closed loops by connecting Wilson 
lines at different transverse coordinates at = ±00. However, this leads to unnecessary com- 
plications because in the light-cone gauge in which we are working, there is nonzero transverse 
field at infinity. 

^ One can make the correspondence explicit by introducing angular coordinates p oc 
(sin^cos^, sin^sin^, cos^). Under this parametrization, p ~ cos^, q ^ (f) and Dp ~ DpDq. 
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with p{x) defined by Eq. (2.15) with J J{x). Correspondingly, Eq. (2.13) is generahzed 
to 

J Dp±^{x)Wr\poo, P-oo] Dp{x+, x) ^ DA^ exp (z-SymI^I + iSwlA', p]) . 

(2.23) 

Eq. (2.23) is the starting point of quantum evolution. By extending the definition of 
gauge transformation to the charge sector, Eq. (2.18), one can explicitly maintain the 
gauge invariance of the theory. Namely, Eq. (2.23) can also be written as 

Z = J Dp±^Wr[poo, P-oo] jj^ Dp{x+) ^ DA'^ exp (iSYuiA^] + tSwiA' , p]) . (2.24) 
2.2 Induced charge from the source term 

In order for Eq. (2.22) to be an acceptable source term, first we have to show that it 
gives the correct induced charge (5p[a~] under one step of quantum evolution A~ — a" + 
SA~. Namely, -S'wfa" + SA~,p] ~ — (p + 5p[a~])5A'. At first sight this looks impossible 
because S'w is linear in A" (so a~ and 6 A" do not couple). However, we observe that the 
corresponding coupling ip'-f^t'^iIjA'^ ~ P"'^a ^^e QCD Lagrangian is also linear in A^. 
As with the perturbative QCD calculation, we expand the exponential e"*''^ in powers 
of A~. To quadratic order, we get 

-^Jdp{x^) J dx+dy+p„(x+)p,(y+)e^^--(M;(x+) + a-(x+))(5A-(y+) + a,-(y+)), 

(2.25) 

Hereafter we suppress the spacial coordinate x = (x^ ,x^). The crucial step of our approach 
is that we rewrite the a~6A~ coupling in Eq. (2.25) as 

- J Dp{x+) J dx+dy+Paix+)pf,iy+)5A-{x+)a;iy+)e'^^^ 
^-j dx+dy+ {e{x+ - y+Wp') + 9{y^ - x+){p'p'^)) 5A-{x+)a^{y+), (2.26) 

where (a = 1,2,3) are non-commutative charge operators. In the zero dimensional 
problem, p" = —gT°'/2. In the three-dimensional case p" are functions of x and satisfy 
local commutation relations 

[p''{x),p\y)] = -^ge-'-p^5^^\x-y). (2.27) 

Eq. (2.26) is the 'magic' of the Wess-Zumino term which can be mathematically justified 
[28]. A path integral of commutative p's endowed with a Wess-Zumino term is equal to 
a matrix element of non-commutative charges p. The ordering of p's follows from the 
ordering of the x"*" coordinate of p{x'^) under the path integral. 
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By decomposing the product of two p's into the symmetric and anti-symmetric parts 

P''p' = \\p\p'] + \{p\p'}, (2.28) 
we write the second hne of Eq. (2.26) as 

(2.29) 

The symmetric term oc {{p"',p^}) in Eq. (2.29) vanishes because it is proportional to 

dy+a^{y+) = a;{p- = 0) = 0. (2.30) 



Eq. (2.30) is vahd since the semihard field is nearly an on-shell excitation [17] having 
A+ > p+ > 6A+ and p\/2bA'^ > p~ > p]_/2K'^, where pj_ is a typical transverse momen- 
tum. [In the LLA, the precise value of p_L does not matter.] Returning to the path integral 
representation {p'^) — * / Dp p'^e^^^'^^ we re-exponentiate Eq. (2.29) and read off (a part 
of) the induced charge from the coefficient of —i5A~ 



£abc „ 

5p«(a;+) ^ -^-^P-oo J dy+{e{x+ - y+) - 9{y^ - x^))al{y^). (2.31) 

The x~^ coordinate of p*^ can be chosen freely, since there is only one factor of p left in 
Eq (2.29). This really does not matter; as we shall argue at the beginning of Section 
3, p is almost static in the JIMWLK regime. In the above, we have set x'^ = — oo. 
Eq. (2.31) agrees with the literature [16, 17], although the derivation here is very different. 
As remarked above, our derivation is closely tied to the corresponding pQCD calculation. 
Indeed, pQCD, the a'SA" coupling arises due to second order perturbation 

_^(^^+ta^M-)(V^7+tVa6-) - -\p''5A-p'a^, (2.32) 

We see that the Wess-Zumino term is a convenient trick to replace the dynamical fermionic 
field with classical charges while correctly retaining the color commutator which in pQCD 
was implied by the quantum commutator. 

Next we consider the virtual contribution to the induced charge coming from the ~ 
pa''a~SA~ coupling. This time we expand up to the cubic order 



^ '^jDp{p{a- + 5A-)fe^ 

IjDpj dx+dy+dz+p'^{x+)p\y+)p%z+)6A-{x+)a;{y+)a:{z+)e''-^^ 



3! 



J dx-^dy^dz-^SA^ (|/"^)«c {z'^)(^xyz{paPhPc) + Oxzy{paPcPb) 

+Oyxz{PbPaPc) + OyzxiPbPcPa) + ^zxyiPcPaPb) + Ozyx{PcPbPa) ), (2.33) 
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where 



e,y, = e{x+ -y+)e{y+ -z+). (2.34) 

The coefficient of —iSA" is another contribution to the induced charge Sp^'^\ For this 
part of the induced charge, all we need in the later developments is the expectation value 
{Sp^"^^) where the averaging (...) is computed in the Gaussian approximation using the 
background field propagator Eq. (3.23). Anticipating this, we perform the replacement 

a;iy+)a:{z+) ^ 5,ciG~iy+ - z+) (2.35) 

already at this point. Note that this component of the propagator is diagonal in color 
indices. This is certainly correct for the free propagator (which is the case for the BFKL 
regime). In the JIMWLK regime, the background field modifies the propagator. Never- 
theless, Eq. (2.35) is valid in practice. See, Appendix C of the third paper in Ref. [17]. 
Since p^Pb is a group Casimir, it commutes with p^. Using this fact and Eq. (2.30), we 
obtain after some algebra 

-^^^ / dx+dy+dz+SA-{x+)iG—{y+ - z+)pa(e{y+ - x+)e{x+ - z+) + e{z+ - x+)e{x+ - y+)). 

(2.36) 

Prom this we can read off the (expectation value of the) induced charge 

(<5p(^)) ^ ^P-ooJ dy+dznG-{y+ - z^){9{y+ - x+)^(x+ - z+) + 0{z+ - x+)9{x+ - y+)), 

(2.37) 

in agreement with the literature. 

Summarizing, we expanded e"*''^" '^^^ ^ and re-exponentiated the relevant (in the sense 
of the LLA) couphngs between a" and SA~. Effectively, we have replaced 



exp^-i J dx^p{x^)(a (x^) + SA (a;+)) 

exp (^-i J dx^p_ooa' - i j dx^{p{x'^) + 5p^^'^ + a'lla"^ , 

(2.38) 

where we again chose p{x'^ — — oo) in the linear term in a~ (see the remark below 
Eq. (2.31)). The last term in Eq. (2.38) is a correction to the background field propagator 
(c.f., Eq. (3.23)) coming from the a'a" term in Eq. (2.25) 

j a-na- ^ ^ / dx^dy+(e{x-^ - y+) - ^(y+ - x+))pi^a„-(x+)a,-(y+). (2.39) 
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3 Quantum evolution and renormalization group at high density 



3. 1 Recovering the CGC picture 



In this section, we revisit the renormahzation group evolution in the high density regime 
[16, 17] with the aim of showing the renormalization of the source term. We start with 
Eq. (2.23) in the hght-cone gauge = and expand around the solution to the classical 
equations of motion which are obtained by varying the field 

D,F^^ = p{x+), D,F^^ = Q. (3.1) 

Eq. (3.1) was solved in [11], with the solution again in the same form as Eq. (2.2) ^ 

A^{x+) = 5'''BHx+), B' = -Ud'UHx+). (3.2) 

9 

Crucial simphfication in the JIMWLK regime is that we can neglect the a;+-dependence 
in Eqs. (3.1) and (3.2). This is because the important paths in the path integral / Dp{x~^) 
are those with weak x'^ dependence. To illustrate this, consider the solution to the saddle 
point equation 

D-p{x^) = 0. (3.3) 
This can be obtained by performing infinitesimal rotation 

S^Se'"', u^uj''^, (3.4) 



in Eq. (2.14) and Eq. (2.15). Under this rotation, p and -Swz change by 

bp" = e^^^uj^p", bSy^z = -- / dx+p"d-u\ (3.5) 

9 J 

Eq. (3.3) immediately follows by imposing — / bp^A'+bST^z — 0. The solution to Eq. (3.3) 
is 

p"{x+) = Wab{x+)pt^ = (iy(x+)p_ooW^^(x+))", (3.6) 

where 

W{x+) = Pexp (^tg J'' dz+A-{z+)^ . (3.7) 



^ In [11], the x+-dependence of p was assumed to be of the form Eq. (3.6). However, the solution 
Eq. (3.2) is more general and is valid for arbitrary a; "'"-dependence of p(x~^). 
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In Appendix A, we derive an interesting property of the saddle point solution Eq. (3.6). 
Since A~ field is weak in the JIMWLK regime, at the saddle point 



^ p_oo. (3.8) 

At high density where p~C(l/5i)^l, the path integral receives important contributions 
from the neighborhood of the saddle point. Therefore, it is legitimate to approximate 
p(x+) p_oo when solving Eq. (3.1). Moreover, 

Wr[poo: P-oo] ~ W^r[P-oo, P-oo] = W^t[P-oo]- (3.9) 

This is a very important approximation which allows us to recover the CGC picture from 
our ^'''-dependent formulation. Namely, as a result of Eq. (3.9), Wr becomes real and 
positive, therefore it literally serves as a weight function. In the high density regime, at 
fixed rapidity, the above argument about the saddle point is generally correct and leads 
to a static, classical theory characterized by the averaging Eq. (2.6). However, the saddle 
point solution and the corresponding approximation Eq. (3.8) have only a limited sense 
when we consider quantum evolution, namely, when the semihard field is introduced. 
Indeed, it was essential to keep track of the a;+-dependence of p{x^) and maintain the full 
(not saddle point) path integral Dp[x^) with the Wess-Zumino term in order to correctly 
derive the induced charge 5p^^''^^[a~\. Therefore, we use the approximation p{x^) ~ p_oo 
only when it is safe to do so. Our criterion of 'safe' is that the semihard field a~ is not 
involved. For example, Wr[Poo,P-oo] is safe but the term —p{x^)a^ {x^) in S'w is not: 
p{x'^)a~{x^) 7^ p-coa~{x^). In this way we can maximally exploit the approximately 
static nature of the charges in the JIMWLK regime without tampering the precise x'^ 
structure of the source term S^ff[p{x'^)\. 



3.2 Renormalization group 

With these caveats in mind, we expand around the static solution B^[p_^] 

[ Dp_^Wr[p-oo] Dp(x+) [ DSA''- 

X exp {iSYM[B' + 5A\ SA'] - i J dx+p6A- + iS^^z[p{x+)]^ ■ (3.10) 

In the above, the field 6A'''^ contains modes with p"*" < A"*". In order to derive the effective 
theory at rapidity t+5t, we decompose the soft field 6A^'~ — > a^'+SA^' and functionally 
integrate over a*'~ 



Z = [ Dp_ooWr[p-oo] r'~^Dp{x+) [ D5A''- r ^ Da''- 
X exp (^iSYM[B' + a' + 5A\ a" + SA'] - i j dx+p{a- + SA') + iSwz[p{x'')]^ (3.11) 
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In the source term we perform the replacement Eq. (2.38). The Yang-Mills action gives 
other contributions [16, 17] to the induced charge 

-5fyYM[a']5A-, (3.12) 

to be added to the induced charge from the source term 

5pYM + 5p^^^ + = 5p[a] . (3.13) 

As was observed in [17], 5p has a support in a narrow strip l/ftA"*" > x~ > l/A'*', while 
the original charges p sit on 1/A+ > x~ > Q (Fig. 2). Therefore, each step of quantum 
evolution piles up a layer of new classical charges at larger positive values of 

The expansion of 

Sym[B' + a' + 5A\ a- + SA'] (3.14) 

requires a care because in the presence of quantum fluctuations the average field is not 
just B'-, but 5* + {5A\^^^[5p]) [17] where 5Ajjj^[5p] is that part of the soft field 5A^ induced 
by 5p and obeys the Yang-Mills equation with the renormalized charge^ 

. =D,F^^'\ . =5^'Hp_^ + 6p\a]). (3.17) 

For given 5p, the solution of Eq. (3.17) is again a pure gauge (see Eq. (2.2)) 

B' + 5A\^^ = ^Ud'U\ W = Pexp (^ig J dx'ia + 5a)^ , (3.18) 

where 5a is the induced field in the Coulomb gauge. Taking this into account, we write 

5A'^5Al^^[5p]+5A\ (3.19) 
and expand Eq. (3.14) around B^ + 5^ijj(j[(5p]. To hnear order in small fields, we get 



^ In the LLA, only the transverse p = i components are important. Indeed, the one- and the 
connected two-point functions of (5^-^^ are logarithmically enhanced [17] just like the one- and 
the two-point functions of 5/3[a]. By expanding the left hand side of Eq. (3.17) to quadratic 
order in ^-^l^^, we obtain relations between the correlation functions: 

" " ^^^['^^^ = I B^^^^^) + I sIIZa^ I B^^K^ - -s5r, (3.15) 

from the p = + component of Eq. (3.17) and 

^ I b(*^».> + kjimjM 1 = 0. (3.16) 

from the p = i components. 
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Sym[B' + SAi^^ + a' + 5A\ a' + 5A-] 
= Sym[B' + MU] + D^F'^+lsi+sAi Xa- + SA-) + D,F^%,^sa^ M + SA') + ■■■ 

md ind 

= (p_oo + Sp[a] ) (5 A- + a-) + • • • , (3.20) 
where we used 

Sym[B' + 5AI^]^0, (3.21) 

because + SAl^^ is a pure gauge. Note that there is a mismatch between Eq. (3.20) and 
Eq. (2.38) (with the total induced charge Eq. (3.13)). The term dp[a]a~ in Eq. (3.20) is 
not cancelled by the source term. This term will play an important role below. 

Collecting all factors, we are now prepared to integrate over the semihard field. In the 
Gaussian approximation, Eq. (3.11) becomes 



Z^l Dp_^Wr[p_^] I Dp{x+)D5A''-Da''- 



X exp(^iSYM[B' + SAi^^[Sp[a]] + SA\ SA-] + ^-a^G-la" 

M5p[a]a~ - i{p + Sp[a])SA~ + iSwz[p] ) , (3.22) 

where is the background field propagator 



We introduce the following trick. 



+ V5,_n. (3.23) 



Z^J Dp^^DSp-^Wr[p-oo] J Dp{x+)DSp{x+)DSA'^-Da''- S{Sp-Sp[a]) 



X expUsYulB' + SAl^[6p] + 6A\ SA^] + -a'^G^Ja" 



+i5pa —i(p + 5p)5A + iS'wz[p] 
= j Dp_^D6p_^WAp_^] j Dp{x+)D6p{x+)DSA''-Da''-D7T 
X expl^iSYM[B' + SAi^^[Sp] + SA\ SA'] + '-a^Q-Ja" 

+idpa~ + i7r{dp - dp[a]) - i{p + Sp)SA~ + iSwz[p]j ■ (3.24) 

The integration over a*'~ can be done. To the order of interest (retaining terms enhanced 
by 6t), we obtain [16] 
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Z = y L>p_oo£>5p_ooW^r[p-oo] / Dp{x-^)DSp{x+)DSA'^-D7T 

-i{p + Sp)SA- + iS^^z\p] + iSwz[Sp\^ , (3.25) 

where a = {Sp[a]) oc agSr and x = (^pN^pH) c< agSr. Note that, in addition to the usual 
building blocks of the JIMWLK kernel, a and x, we have included a Wess-Zumino term 
for the induced charge Sy/z[Sp{x^)]- This term does not contain a logarithm, but must 
be included in order to render the theory at rapidity t + 6t gauge invariant. [The sum 
—5p5A''+S\^z[Sp] is gauge invariant.] We are not able to present a complete derivation of 
S'wz [5p] here, but argue how this term could arise. The origin of S'wz [^p] is the elusive Spa~ 
term. However, it does not arise in the Gaussian approximation. Indeed, if we performed 
a Gaussian integral in Eq. (3.24), we would get an insignificant term 

1 dxdy6p\x)G-,-{x,y)5p\y), (3.26) 

which is not the Wess-Zumino term, [Remember that, despite our notation, S^ffz[5p] 
cannot be written explicitly as a functional of 5p !] nor does it contain a logarithmic 
enhancement. ^ The absence of S'wz[^p] in the Gaussian approximation is understandable. 
As seen in Eq. (2.14), the Wess-Zumino term is defined over the entire gauge group. 
Therefore, its origin must have to do with gauge invariance. Once we truncate the Yang- 
Mills action in the Gaussian approximation, invariance under large gauge transformations 
is lost. This suggests that in Eq. (3.24) we have to perform an integral with the full Yang- 
Mills action 

f[5p] = j Da^' e:!L-^{iSYu[a\cr] + i j dx^5pa'^, (3.27) 

where we neglected the external fields and p. By this we formally treat the fields a*'~ as 
large and assume that the external fields arc not essential for the generation of S^z{^p\- 
[S'wz[^p] does not depend on or p.] Indeed, it is known [30] that the Wess-Zumino 
term comes from the large a*'~ region of the integral Eq. (3.27). To see this, consider the 
following change of integration variables in Eq. (3.27) [See [30] for more details.] 

a'- a'- = V^a'-V + V^9^'-y", (3.28) 

where V does not depend on x~ . This is a gauge transformation which preserves the 
light-cone gauge condition a+ = 0. Using the gauge invariance of the Yang-Mills action, 

^ Note that the x'^ dependence of 5p and G originates from that of the semihard field a^. 
Therefore, the p" integral {p\ / 26A+ > p" > p]_/ 2A+ ) in the Fourier representation of G (x, y) 

contains an oscillating phase e^^^ ~^ which prohibits a logarithm. This is in contrast to 
Eq. (2.37) where the integration over y+ and z'^ can be explicitly performed and a logarithm 
does arise. 
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one can easily check that 



f[5p] = f[V5pV^]e-^P^^^~^. 



(3.29) 



A general solution for Eq, 



(3.29) is 



f[5p] = ^[5p]e''^wz[M 



(3.30) 



where g[5i^ is a gauge invariant function 



g[5p]=g[V5pV^]. 



(3.31) 



In Eq. (3.30), the Wess-Zumino term arises as a special solution to Eq. (3.29). The rele- 
vance of this term is easy to understand. The extra phase factor in Eq. (3.29) due to the 
non-gauge invariance of / 5pa~ is precisely generated by the variation of the Wess-Zumino 
term. The function g[5p] cannot be calculated exactly, but since the Gaussian approxima- 
tion Eq. (4.2) (with G the free or dressed propagator) is already unimportant (i.e., not 
enhanced by a logarithm), we can simply set g = 1 in the leading logarithmic approxima- 
tion. This is the reasoning of our addition of the Wess-Zumino term in Eq. (3.25). 

Returning to Eq. (3.25), we observe that 



This holds because p and 5p have different supports in x [17]: 1/A+ > x > for p, and 
l/foA"*" > x~ > 1/A+ for 5 p. Moreover, from Eq. (3.25), we can easily deduce the evolution 

of Wr (p'_oo = P-oo + (^P-oo) 



in agreement with the path integral formula previously derived in [31]. ^ ^ The infinitesi- 
mal evolution Eq. (3.33) is equivalent to the JIMWLK equation 



where the subscript r in p means that the derivatives are taken at the highest value of 
x^] = 1/A+ [17]. [Remember that the support of 5p_oo in Eq. (3.33) is l/hlv^ > x^ > 
l/A"*".] Using Eq. (3.32) and Eq. (3.33), we finally arrive at 

^ Here we neglect the x"*" dependence of 5 p. See the remarks following Eq. (3.9). Eq. (3.33) 
shows that 5p_oo is Gaussian distributed with the mean a and the variance x (which are both 
static). 

^ See, Eq. (4.11) of Ref. [31]. There the authors work in the Coulomb gauge Eq. (2.4). There is 

a subtlety in going from the light-cone gauge to the Coulomb gauge. The gauge function which 
realizes this rotation is not Eq. (2.3), but Eq. (3.18). See the discussion in Section 3 of the first 
paper in Ref. [17]. 



'S'wz[p] + S^zih] = S^z[p + M- 



(3.32) 




(3.33) 




(3.34) 
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X exp [iSYM[B" + 6A\ M"] - ip'M" + iSwz[p' {x+)]) , (3.35) 

where B'^ — B'^ + SAl^^[Sp]. Eq. (3.35) is exactly the same form as our starting formula 
Eq. (3.10). This is the JIMWLK renormahzation group picture. Namely, the quantum 
effects at one step of evolution is completely absorbed by the change of the weight function 
and the form of the effective theory is preserved. 



4 Conclusions 

In this paper, we have investigated the role of the Wess-Zumino term in the high energy 
limit of QCD where the target becomes a Color Glass Condensate. We have shown that 
the Wess-Zumino term can be naturally incorporated in the JIMWLK formalism as a 
part of the source term Sy^. The simple eikonal coupling between the gauge field and the 
source —p5A~ makes the renormahzation group description more transparent. We have 
argued that after one step of quantum evolution the functional integral Eq. (3.10) remains 
the same form including the source term. This last point has not been discussed in the 
existing proofs of renormahzation group [16, 17]. 

Finally, we expect that the real strength of our approach is that it can be straightforwardly 
applied to other limits of high energy QCD, namely, the dilute regime (Fig. 1(b)) and the 
Pomeron loop regime (Fig. 1(c)), where the importance of the Wess-Zumino term has 
been first recognized. The color glass averaging in the dilute regime proposed in [6] is (see 
also [8]) 

(pVp^ • • • )r = / Dp{t)WMt)Mti)p\t2)p%h) (4.1) 

where the ordering variables ti are constrained such that ti > t2 > ^3 > " " " , but otherwise 
arbitrary. The "weight function" VrT[p(^)] contains a complex phase e''^wz[p(*)] which en- 
sures the color commutators of p. It follows from our analysis in Section 2 (see Eq. (2.23)) 
that Eq. (4.1) should read 

•••).= / Dp±^Wr[p±oo] r L'p(x+)e^^wz[p(-+)]p«(^+)pfe(4)pc(^+) . . . ^ (4 2) 

with the same Wess-Zumino term appearing in Eq. (2.22).^*^ Therefore, the JIMWLK 
formalism with our new source term allows us to treat different limits of high energy 
QCD in a single framework. The renormahzation group description of quantum evolution 
in the dilute regime is subtler than in the JIMWLK regime because the division between 

Note that the approximation Eq. (3.9) is invalid in the dilute regime. This means that VF,-[p±oo] 
is complex even without the Wess-Zumino term. 
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the "quantum" and "classical" theories is shifted towards the quantum side. [Compare 
Eq. (4.2) with Eq. (2.6).] We leave this problem for future works. 
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A Remarks on the saddle point solution Eq. (3.6) 



For possible future applications, in this appendix we note a curious property of the saddle 
point solution Eq. (3.6). First we observe that the source term vanishes at the saddle 
point 



= 0. (A.l) 

saddle point 



To show this, return to Eq. (2.15) 



We see that S{x~^) and S'_oo are related by 

S{x+)^S_^W^{x+). (A.3) 
Therefore, the Wess-Zumino term Eq. (2.14) becomes 

Swz[p{x+)] .=^j[ dx+tr [r^S^^W^'d-WSlj 

saddle point J 

= -gj j dx^trlT^S^^W^T^WSi^]^ = J dx+ p''{x+)5A- . (A.4) 

Since S'w is gauge invariant, Eq. (A.l) holds in any gauge. Then let us consider the 
exponential factors in Eq. (3.10) in the Coulomb gauge Eq. (2.4) 

SYM[a, SA-] - I pSA- + 5w[p(x+)], (A.5) 
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where we neglected SA^. In the context of [11], Eq. (A. 5) should be called the tree level 
effective action. Since the only non-vanishing components of the field strength are 
the Yang-Mills action becomes 

_1 j pi^^F^, ^ j F+'P-' = j d'ad'SA- ^- j {V''a)5A- = (A.6) 

Therefore, in the Coulomb gauge, at the saddle point Eq. (3.6), all the three terms in 
Eq. (A. 5) are equal up to a sign. In particular, 

-S'ym = 'S'wz- i^-"^) 

Eq. (A. 7) was conjectured in [6] as a possible origin of the Wess-Zumino term. Although 
at the moment we do not see a connection between our derivation and the argument in 
[6], in any case it is interesting to pursue physical implications of Eq. (A. 7). 
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